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Abstract 

We study the properties of a thermodynamic system having the 
symmetry of a quantum group and interacting with a harmonic po- 
tential. We calculate the dependence of the chemical potential, heat 
capacity and spatial distribution of the gas on the quantum group 
parameter q and the number of spatial dimensions D. In addition, 
we consider a fourth-order interaction in the quantum group fields ^, 
and calculate the ground state energy up to first order. 
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1 Introduction 



The recent observations of Bose- Einstein condensation [|T], ^ has triggered 
numerous theoretical investigations on the behavior of a trapped Bose gas 
under experimental conditions. Furthermore, the stable fermionic isotope of 
lithium, ^Li, has been trapped and cooled in a similar way as its bosonic 
counterpart ^Li [||. One of the most interesting aspects of a trapped Fermi 
gas is that, as reported in References p, for a sufficiently strong trap- 
ping magnetic field a BCS transition is expected to occur at densities and 
temperatures experimentally attainable. The critical temperature and one- 
particle excitations for the superfiuid phase transition have also been studied 
[0, ||. Two recent articles |T0| on trapped Fermi gases, in which the BCS 
transition is excluded, include a study on their thermodynamic behavior in 
harmonic traps. 

The purpose of this paper is to study the behavior of a thermodynamic 
system having the symmetry of a quantum group. The role of quantum 



groups in physics has its origin in the quantum inverse scattering method [|ll 



and the possible relevance they may play in other fields of physics has been 
the focus of great interest ever since. A quantum group, as compared with a 
classical group, contains an additional parameter q. In the quantum inverse 
scattering method the parameter q acquires a physical meaning through its 
relation with the Planck's constant. In two recent articles [0, |13[ we con- 
sidered the system represented by the hamiltonian TIb = 
2 = 1,2, with the operators and the adjoint $j satisfying algebraic rela- 
tions which are covariant under the quantum group SUq{2). For q = 1, the 
operators and become standard creation and annihilation boson oper- 
ators respectively. In this model, Bose-Einstein condensation is realized as a 
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second order phase transition, and the heat capacity becomes discontinuous 
at the critical temperature in one, two and three dimensions, even without 
the presence of an external potential. In a similar fashion, one can define 
a set of operators \E'i and \E'j, with corresponding quantum group covari- 
ant algebraic relations, such that for g = 1 they become standard fermionic 
operators. 

A quantum group invariant hamiltonian is, at the moment, a mathemat- 
ical model with phenomenological implications that depart from those pre- 
dicted by either boson or fermion gases. This departure is measured by the 
value of the parameter q. There is, at present, no phenomenological evidence 
that quantum group symmetries are realized in any particular thermody- 
namic system. Since the operators ^ satisfy the Pauli exclusion principle 
and their representation in terms of ordinary fermions ip is not linear, any 
hamiltonian written in terms of them can be interpreted as an interacting 
fermion system. Our main motivation focuses in investigating about the role 
that these interactions, required by quantum group invariance, may play at 
low temperatures. The interaction term is fixed by the quantum group and 
its strength is parameterized by q. Naturally, for g = 1, these interactions will 
vanish and our results will reduce to the case of a gas with two fermion species 
interacting with a harmonic potential. Within the context of our results, only 
experiments with trapped Fermi gases can tell us whether quantum group 
symmetry is realized in the behavior of real gases. A comparison with the 
results reported in this paper will indicate whether a particular, and presum- 
ably small, deviation of the parameter q from the standard value q = 1 fits 
the experimental data. Therefore, the calculations considered here, will help 
to answer questions regarding the possible relevance of quantum groups in 
a thermodynamic system. In this paper we study the thermodynamics of a 
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system composed of quantum group fermions interacting with a harmonic 
potential For matter of simphcity, we will consider the case of SUq{2). 

This paper is organized as follows. Section ^ is a very brief introduction to 
the formalism we will use in the following sections. A detailed discussion on 
quantum groups can be found in References flSl . In Section |^ we calculate the 
properties of this model interacting with an isotropic harmonic oscillator. In 
particular, we calculate the chemical potential, internal energy, heat capacity, 
and the particle distribution at T = for an arbitrary number of spatial 
dimensions D. In Section ^ we introduce a fourth order \E'-interaction and 
calculate the ground state energy up to first order. In Section |^ we summarize 
our results and compare them with the g = 1 case. 



2 Free SUq{2) fermion gas 



The quantum group GLq{2, C) consists of the set of matrices T 



with noncommuting elements generating the algebra Aq [16 



a b 
c d 



ah = q~^ba , ac = q~^ca 
be = cb , dc = qcd 
db = qbd , da — ad = {q — q^^)bc, 



(1) 



with the quantum determinant 



detqT = ad — q ^bc. 



(2) 



The quantum determinant is defined such that it belongs to the center of the 
algebra. The relations in Equation (|I]) have the remarkable property that if 
the elements a, b, c,d of T commute with the elements a', b', c', d' of T', then 
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the elements of the matrix TT' also generate Aq. By setting detqT = 1 and 
the unitary conditions [jT^ 

a = d , b = —q^^c, (3) 

leads to the quantum group SUq(2) with q E R. It is clear that the quantum 
group SUq{2) — > SU{2) as the parameter g — 1. In Reference we 
introduced a set of operators that transform under SUq{2) and become 
ordinary fermions ipi in the q —>■ 1 limit. This is accomplished by the set of 
relations 

{^2,^2} = 1, 

{^1,^1} = 1 - (1 - g-2)^2*2, 

vl/^v[/2 = -g\l/2\l>i, 
^^v[/2 = -q\i/2^i, 
{^1,^1}= ={^2,^2}. (4) 

Equations are covariant under the action of the group SUq{2) 

^. = E^..^., (5) 
i=i 

and they become the SU{2) covariant fermion algebra for q = 1. According 
to Equations (P, the operators have a representation in terms of ordinary 
fermions ipi as follows 

^1 = ^1 [1 + (q-^ - l)4i^2\ , ^2 = ^2, (6) 

with the corresponding equations for the adjoint Therefore, the simplest 
quantum group invariant hamiltonian Ti. = J2k ^k^k is written in terms of 
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fermion operators as follows 

rL^Y.En [^in^l,n + V'LV'2,« + - 1)V'LV'1aV'2,«V'2,«] (7) 

The grand partition function is given by 

K 

which for g = 1 becomes the partition function for a gas with two fermion 
species. The occupation number (n) is a function of the energy according to 

HE)) = ^ (9) 

where z is the fugacity and the function 

For q > 1 the function {n{E)) is very similar to the Fermi function. The 
chemical potential at zero temperature is independent of the number of space 
dimensions D 

Mq > 1) = ^-^/^o, (10) 

where fio is the Fermi energy for q = 1. For q < 1, {n{E)) departs consider- 
ably from the Fermi case. At T = 0, states with energies < E < q^^Q{q) 
are fully occupied, and those with energies q'^fJ^oiq) E < fio{q) have occu- 
pation numbers (n) = 1/2. The chemical potential, at T = and g < 1, is 
related to the Fermi energy for = 1 as follows 

/ 2 Y^"" 

l^oiq < 1) = J ^^^^ 

The total number of particles, in the thermodynamic limit, is given by the 
integral 

roo 

(TV) = 2/ g{E){n{E))dE, (12) 
where g{E) is the density of states. 



3 SUq(2) fermion gas in a harmonic oscillator 

The single-particle levels, excluding the zero-point energy, are given by the 
well known formula 

D 

E = huY,nj. (13) 

i=i 

As an approximation, , we consider a continuous spectrum with density of 
states 

A more accurate expression for the density of states was provided in Refer- 
ences. [0 and applied to the Bose-Einstein condensation of a finite number 
of bosons in a harmonic oscillator trap, and in Reference [|ry] for the study of 
a free fermion gas interacting with a harmonic oscillator potential. Equation 
([mD will be sufficient for our purposes. 



3.1 Chemical potential 

The chemical potential, for low temperatures, can be obtained by splitting 
the integral in Equation (|I2[) in the following intervals 



a) g< 1: [0,q^fi{q)], [q^M,fi{q)], [^(g),oo], 

b) g> 1: [0,2/x(g)/(l + g-2)], [2fi{q)/{l + q-^,q^fi{q)] , [gV(g),oo] 
such that reverting the solutions leads to the results 
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l + q' 



[D-l] 



1.64 

+ 



-2^2 



'kT\ 



(15) 



A simple inspection shows that for q < 1 the chemical potential depends 
also linearly on the temperature. In contrast to the fermion case, where 
/i is constant for D = \, the function yu(g < 1) decreases with T. Figure 
1 shows a graph of the chemical potential for q = 0.2, 1, 2, obtained from 
a numerical calculation of Equation (Jji^. The chemical potential at zero 
temperature, ^o{q), is independent of D for g > 1. For q < 1, /io(g) decreases 
as D increases. For low temperatures, the function fi{q < 1) depends more 
strongly on q as we reduce the number of dimensions D. 



3.2 Internal energy and heat capacity 

The internal energy is obtained from the partition function 

d\nZ 



U 



dp 



+ /i(iV) 



f[z,q) 



(16) 



The integral in Equation (16) can be solved for low temperatures by using 
the same kind of approximation we used to obtain the chemical potential. 
For g 7^ 1 the internal energy is given by the equations 



U{q < 1) 
(AT) 



D 



2(^+1)/ 



D 



2^"(D + l)(l + g2i3)i/i3 



+ f /^o (2.88 - |i _^ 



2D\2 \ 

^ ^ In^S 



2DU/D /f^r^ 



21/D 



— (17) 



U{q > 1) 
(AT) 



D 



l + q- 
D + 



1.64 (kT" 



l + q- 



(18) 
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Figure (2) shows a graph of the heat capacity 

/ ^ -.U ^dU f da u\ , . 

= (Z, + l)- + .,3_(^^-^j (19) 

for D = 3 and q = 0.2, 1,2 which results from a numerical calculation of 
Equation (]TP|). 

3.3 Particle distribution at T = 

In the semiclassical approximation, Thomas-Fermi approximation, the spa- 
tial distribution is given by the integral 

"'-'■■^^'nmi^^ (20) 

Replacing Equation and defining the variable K = p^/2m, this integral 
becomes 

^ ' ' T{D/2) {27rnf Jo i + 2e-/3(^-A) + e^+A-a+^-^W 

(21) 

where for convenience we have defined 

A = - ^muj^r'^, (22) 



and 



fx = /i(g) —mtu'^r'^. (23) 



For g > 1 we have /io > jxo, and Equation (pT]) for T = can be easily 



calculated if we divide the integral in the following intervals: 

[0, /io], [Ao, (/to + /io)(l + g"^)], [(Ao + /io)(l + q''^),q^h] and [g^/io, oo] 
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leading to the result 



L>/2 



n(r,0) = 2A 1 - 



, g > 1, 



(24) 



where 



A 



27r^/^ {D-1)\{N) 
T{D/2) 27iDr^ ' 



and rp = y2/io/mc(j^ is the Fermi radius for g = 1. This equation shows 
that the spatial distribution is independent of the parameter q for all the 
values q > 1. A numerical calculation of Equation (^T]) shows that this 
independence on q remains valid for T > 0. 

For g < 1, a similar calculation leads to the result 

D/2 




2D\ ^/D^ 



1 



l + q 



2D\ ^/D^ 



, /n\ D/2 



(25) 



(26) 



for r > 



qrp- From this results we see that the gas is 
less confined for g < 1. This fact becomes an expected result after rewrit- 
ing the hamiltonian in terms of fermion operators. The fermion interaction 
term induced by the requirement of quantum group invariance becomes more 
repulsive as we decrease the value of q below q = I. 



D/2' 
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4 Higher order interactions 



In the previous section we studied the simplest quantum group invariant 
system interacting with a harmonic potential. According to equations 
introducing interactions between the operators will lead to higher order 
interactions terms when the hamiltonian is rewritten in terms of ordinary 
fermion operators. We consider the simple interaction hamiltonian 

Hi = g ^p;,i^P^,2^p2,2^pi,i, (27) 

where {"^p'^i, ^pj} = 0, Vz, j and p' ^ p, and quantum group matrix elements 
{a{pi) , b{pi) , c{pi) , d{pi)) commute with the set {a{pj),b{pj),c{pj),d{pj)) for 
i ^ j. In equation (|27| ) the coupling g is related to the scattering length a 
by 

and the summation is over all momenta satisfying 

Pi + P2 = Pi + P2- (29) 

The first-order correction to the energy of the system is given by those terms 
with p'j^ = pi and P2 = p2. Again, it is convenient to rewrite the hamiltonian 
in terms of ordinary fermions. Replacing Equation (|^) in Equation ( P?] ) leads 
to the first-order correction for the energy 

E^'^ =9^-^+ 9{<1-' - E ^..1^.-2, (30) 

where Np^i = i/jji/ji and we took the equilibrium values 

EA^P. = ^- (31) 
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In the ground state, all the fermions occupy the states with lower momenta, 
such that the sum J2pNp,iNp,2 = {N)/2. Thus, for the ground state first 
order correction we obtain the simple expression 

4"=™-^^. (32) 

Equation (|3^) shows that the first order correction to the ground state energy 
decreases as we increase the value of the parameter q. From Equations (p!7|) 



and ([T8|) for T = we write, up to first order, for the ground state energy 

^o(^>l) = ^IL/n + + (34) 

Figure 4 shows the ground state function Eq = 4:Eo/g{N)^ for D = 3. We 
have assumed the case of a dilute gas, na^ = 10"^, with scattering length 
a ^ 5 X 10~''cm, trap size L ^ 5 x 10~^cm and m = lO'^'^^kg. The value 
of the parameter q is restricted, for g < 1, by the condition Eq^^ < Eq^\ 
For the case displayed in Figure 4 the range of allowed values is given by 
0.09 < g < oo. 



5 Conclusions 

In this paper, we first studied the thermodynamics of the simplest quantum 
group invariant hamiltonian including a harmonic potential. For q = 1, this 
system becomes a fermion gas with two degrees of freedom in a harmonic 
potential trap in a D-dimensional space. We have analyzed this model by 
using the representation of the operators \E'j in terms of fermion operators. 
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A calculation of the chemical potential shows that for g < 1 it acquires a 
linear temperature dependent term no present in the fermion, SU{2), case. 
One consequence of having this linear term is that for D — 1 the chemical 
potential is not constant but decreases with the temperature. At moderate 
temperatures the heat capacity increases with q, and at high temperatures 
C = Dk{N), becoming independently of the value of q. A calculation in 
the Thomas-Fermi approximation at T = shows that the particle spatial 
distribution is independent of g for g > 1. For g < 1, if we decrease the 
value of g, it makes the interaction more repulsive and thus the gas becomes 
less confined. Since this fermion interaction vanishes for g = 1, the depletion 
of the gas for g < 1 is a direct consequence of imposing quantum group 
invariance. In the last Section we included a fourth-order interaction in 
terms of the ^ operators, and calculated the ground state energy up to first 
order. This fourth-order interaction becomes a sixth-order interaction in 
terms of fermion operators. A graph of the function Eq illustrates that even 
small deviations from the standard value q — 1 have a nontrivial effect on 
the ground state energy. The ground state energy decreases as the value of 
g increases and it becomes approximately constant for g > 5. 

Our main goal has been in studying the phenomenological aspects of a 
model having the quantum group SUq{2) as its symmetry group. In the same 
fashion that a quantum group is considered a deformation of a classical Lie 
group, we could consider the SUq{2) symmetry of our model as the result of 
a small breaking of SU{2) spin symmetry. An alternative view would be to 
regard SUg{2) as an additional symmetry independent of the spin degrees of 
freedom. 
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FIG. 1. The chemical potential ^{q) for D = 3 as a function of T for the 
parameter values q = 0.2, 1, 2. 

FIG. 2. The heat capacity, as given in Equation (|19|), for D = 3 and the 
values q = 0.2, 1, 2. 

FIG. 3. The spatial distribution for = 3 as a function of r/rp- 

FIG. 4. The dependence of the function = 4:Eo/ g{N)'^ on the parameter 
q for na'^ = 10~^, scattering length a ~ 5 x lO^^cm and trap size L ^ 
5 X 10~'^cm. 
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